Single nonmagnetic impurity resonance in FeSe-based 122-type superconductors 

as a probe for pairing symmetry 
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The effect of a single nonmagnetic impurity in A !/ Fe2-iSe2 (A=K, Rb, or Cs) superconductors 
has been studied based on a three-orbital model. The local density of states on and near the 
impurity site has been calculated by solving the Bogoliubov-de Gennes equations self-consistently. 
Both repulsive and attractive impurity scattering potential are considered in our calculations. The 
impurity-induced in-gap bound states are found only for attractive scattering potential in the state of 
d x 2_ y 2 wave pairing, and it turns out that they are very sensitive to the magnitude of the scattering 
potential. The emergence of the impurity-induced bound states in the vicinity of the Fermi level 
demonstrates a strong violation of the electron-hole symmetry which is originated from the nodeless 
d x 2_ y 2 wave pairing state. The results obtained in our calculation, which simulate the doping of Co 
and Ni in FeSe-based 122-type superconductors, as an approach to examine the pairing symmetry 
of this novel superconducting material, can be a proposal of STM observation. 

PACS numbers: 74.70.Xa, 74.25. Jb, 74.20.Pq 



I. INTRODUCTION 



The FeSe-based 122-type A y Fe 2 - x Se 2 (A=K, Rb, 
or Cs) high-Tc superconductors have been focused 
intensively due to its distinctive high transition 
temperatures(~30K)i^ and novel Fermi surface 
topology. Angle-resolved photoemission spectroscopy 
(ARPES^ and density functional theory (DFT) 
calculations 1 ^ - — show that there is no hole pocket 
around the T point of the Brillouin zone (BZ), and 
large electron pockets are observed around the M point 
which is located at the corners of the folded BZ S . Effort 
has been made to investigate the pairing symmetry of 
the Fe-based superconductors. The popular s± pairing 
symmetry in iron pnictide superconductors in which 
the hole pockets are at the T point, has been predicted 
by spin-fluctuation theor y 14 ' 15 . However, for FeSe- 
based 122-type superconductors, due to the absence of 
the hole pockets, nodeless superconducting (SC) gap 
can be originated from various pairing symmetries. 
Spin-fluctuation theory predicted d x 2_ y 2 wave pairing 
stated— , while s wave pairing state was argued to be 
more favorable against strong disorder—. Moreover, a 
co-existence candidate of anisotropic s and d x 2_ y 2 wave 
pairing symmetries has also been predicted^ 1 -. These 
results reveal the competitive nature of the pairing state 
in these novel superconductors. 

The response of a SC system to the impurities can be 
a probe for the SC pairing symmetry —. The local elec- 
tronic states can be influenced by impurities evidently. 
In a multi-band SC system, the SC order parameter may 
or may not be sensitive to the existence of impurity 23 . 
Doping of Co and Zn in iron pnictide superconductors has 
been realized by several experimental groups 2 ^—. Theo- 
retically, the effects of impurity in iron pnictide supercon- 
ductors has been investigated intensively 2 ^—. The cal- 
culated results of T-matrix method based on a four-band 



model shows the in-gap impurity resonance states^ for 
s± wave pairing symmetry, and results obtained by solv- 
ing the Bogoliubov-de Gennes (BdG) equations in the 
framework of a detailed five-orbital model also demon- 
strates that the s± wave pairing symmetry is character- 
ized by a bound state in the vicinity of the Fermi level in 
the presence of a single nonmagnetic impurity 34 . How- 
ever, for FeSe-based 122-type superconductors, the impu- 
rity resonance bound state, which is located away from 
the Fermi level, were obtained only in d x 2_ y 2 wave pair- 
ing stated. 

In this work, we study the effect of a single nonmag- 
netic impurity in FeSe-based 122-type superconductors 
based on a three-orbital model and by using BdG 40 equa- 
tions. The local density of states(LDOS) and the real- 
space SC order parameters are calculated. It has been re- 
ported that the effective on-site impurity potentials can 
be repulsive or attractive which is related to the energy 
difference between the impurity 3d levels and Fe 3d level 
for various impurity atoms. Realistic estimation of the 
impurity potential in iron pnictide superconductors pre- 
dicted that the potential difference of Zn, Co, and Ni are 
negative, whereas Mn are positive 3 ^. Therefore both re- 
pulsive(positive) and attractive(negative) scattering po- 
tentials have been considered in our calculation which 
may simulate different dopants, such as Zn, Co, Ni, and 
Mn etc.. Responses of isotropic s, anisotropic s, and 
d x 2_ y 2 wave pairing states to impurity are examined, re- 
spectively. We focus on the difference of the resonance 
states induced by impurity among different pairing sym- 
metries and the properties of the impurity-induced bound 
states. 

This paper is organized as follows. Section II presents 
the Hamiltonian of the three-orbital model and the self- 
consistent BdG approach. In section III, we discuss and 
compare the results of the properties of the impurity 
states for different pairing symmetries. Finally, a sum- 
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mary is given in section IV. 



II. MODEL AND METHODOLOGY 

The Hamiltonian of the impure SC system is of the 
following form 

H = Hq + Hp a i r + Hi mp (1) 

and the hopping, pairing, and impurity terms read 

H pa i r = ( A ?f a Lt a ]pi + h - c -) 

Himp =U S a laa a lof (2) 
a, a 

where o\ a!y {ai ai y) denotes the creation(annihilation) op- 
erator of electrons with spin a =7, \. and orbital a at 
site i. t?j are the hopping integrals and /i is the chem- 
ical potential. A single nonmagnetic impurity is located 
at site I. Results of DFT calculation demonstrates that 
the electronic structure in the vicinity of the Fermi level 
is dominated by d xz , d yz , and d xy orbitals 13 , there- 
fore it is feasible to describe the band structure by in- 
troducing an effective three-orbital model^. The hop- 
ping parameters which are given in irreducible hopping 
subsets^ are shown in the Appendix. Fig. Q] shows 
the heavily electron-doped Fermi surface obtained with 
a chemical potential /i=0.312 eV corresponding to a fill- 
ing factor n=4.23. The SC gap function stemming from 
the mean-field decoupling of the pair scattering process 
ytj 13 a ia\ a jPi a 3Pi a ia^ is expressed as 

A ff = V^(a j0i a lat ) (3) 

Single static impurity in a uniform superconductor 
which manifests itself as a point-like imperfection can be 
regarded as a localized orbital-independent on-site scat- 
tering potential U s S^ iRg with the impurity at Rq in or- 
bital representation, provided that the Coulomb interac- 
tion is screened at length scales comparable to the lattice 
spacing 22 . Both positive and negative scattering poten- 
tial U s are considered in our calculation, and by checking 
the results obtained from several different values, it has 
been found that U s = ±6.0 eV can be regarded as the 
unitary limit. 

The Hamiltonian of the impure system which is orig- 
inally written in orbital representation is of its diago- 
nal form in qusiparticle representation by conducting the 
Bogoliubov-Valatin transformatio n 40 ' 41 
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FIG. 1: (color online) The Fermi surface in the extended BZ. 
The nodal lines of the gap function for anisotropic s and 
d x 2_ y 2 wave pairing symmetry are plotted with green and 
red lines, respectively. 

where the quasiparticle creation operator 7^ CT is the lad- 
der operator of the eigen-spectrum of the Hamiltonian 
which satisfies {y\ a ,H]- = -e„7L- 

The diagonal condition of the Hamiltonian is the BdG 
equation 

S(##)(*M4) ,5) 

where 

hff = t"f - ^SijSajs + U s 5uSjiS a p 
and the gap equation is represented as 

n 

The gap function of BCS theory 4 ^ is an eigenvalue 
equation of pairing potential in k-space^, the symmet- 
ric property of which is characterized by the irreducible 
representation of the symmetry group of the Hamilto- 
nian. Therefore, the symmetry of the gap function is 
determined by the form of the pairing potential. The 
absence of the hole pockets of the FeSe-based 122-type 
Fermi surface around T point eliminates the possibility of 
nodes. Calculation pertaining to the magnetic exchange 
couplings reveals that the leading pairing is originated 
from the intra-orbital pairing, whereas the inter-orbital 
components are found to be significantly small^. Con- 
sequently, the s and d x 2_ y 2 wave pairing symmetries are 
considered for intra-orbital pairing only in order to inves- 
tigate their responses to a single nonmagnetic impurity 
scattering. For isotropic s wave pairing, 

Vg* = -.90% (7) 
for anisotropic s wave pairing, 

v ij — ~-^\°i-\-x+y.j ~r Oi-x+yj 

H~ ^i — x — y.j ^~ ^i-\-x—y,j) (p) 
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and for d x 2_ y 2 wave pairing, 

+ 5i- £!j + 5i-$,j) (9) 

where 30,1,2 are pairing amplitudes for each pairing sym- 
metry and x(y) denote the unit vectors along x(y) di- 
rections of the lattice coordinate. The pairing pat- 
terns in real space result in specific pairing symmetries 
of the gap function in k-space as for isotropic s wave 
channel A Q "(fc) ~ A , for anisotropic s wave channel 
A aa (k) ~ cos(k x ) cos(ky) , and for d x 2_ y 2 wave channel 
A aa (k) ~ cos(k x ) — cos(fcy) in the extended BZ, respec- 
tively. The pairing amplitudes are taken to be go = 0.3, 
gi = 1.0, and go = 0.5 which guarantees that the SC 
coherence length £ < 4a, where a is the lattice constant. 

The energy spectrum of the quasiparticle i.e., the 
LDOS at site i is calculated via 

= at- E[KI 2<5 ( e " e «) + K«l 2 ^ e + e «)] ( 10 ) 

J ' orb 

where the prefactor comes from taking the average of 
three orbitals. The Lorentzian smearing method is used 
to visualize the LDOS with a broadening width a = 
0.001. All the self-consistent calculations are performed 
on a 20 x 20 lattice with a periodic boundary condition at 
temperature T= 0.1K, and the single impurity is located 
at the center of the lattice. 



III. RESULTS AND DISCUSSION 

The calculated LDOS for isotropic s wave pairing state 
are shown in Fig. [5]for U s = 1.5 eV and U s = —2.0 eV, re- 
spectively. Obviously, the existence of impurity does not 
change the global properties of the quasiparticle energy 
spectrum. Similar results are also obtained in Ref i 29 i 34 
which demonstrates that there is no in-gap bound state 
for such a pairing symmetry, whereas the order parame- 
ter on the impurity site are suppressed evidently. 

Results of anisotropic s pairing symmetry are shown 
in Fig. |3] for repulsive scattering potentials and in Fig. 
2] for attractive cases, respectively. The SC coherence 
peak is located at ±0.044 eV, and the typical resonating 
states induced by the impurity appear for both cases at 
±0.028 eV for U s = 0.5 eV and ±0.025 eV for U s = -1.0 
eV, respectively. The resonance peaks move towards 
the gap-edge as the amplitude of the scattering poten- 
tial increases for repulsive interaction, while for attrac- 
tive cases, the resonating peaks maintain at the gap-edge 
for small value of scattering potential U s — —0.5 eV, as 
shown in Fig. 0|a) . Additionally, it has been noticed that 
the resonating states for positive interaction U s = 0.5 eV 
is very robust in comparison with other values of U s . At 
unitary limit, the characteristics of the spectrum for both 
repulsive and attractive interaction tends to be identical 
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(a) (b) 

FIG. 2: (color online) The LDOS for isotropic s wave pairing 
state for U s = 1.5 eV (a) and U s = -2.0 eV (b). LDOS for 
impurity, NN, and NNN sites are plotted in red, green, and 
blue lines, respectively. The spectrum of the clean system is 
also shown in black line for comparison. 
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(c) (d) 

FIG. 3: (color online) The LDOS for anisotropic s wave pair- 
ing state for repulsive scattering potential U 3 = 0.5 eV (a), 
U e = 0.8 eV (b), Us = 2.0 eV (c), and U„ = 6.0 eV (d). 



since the electrons cannot hop to the impurity site, re- 
gardless of the detailed feature of the scattering interac- 
tion. We note that the resonating states within the SC 
gap, as compared with the clean cases, appear always in 
pairs which are located symmetrically with respect to the 
Fermi level. The magnitudes of the LDOS of the resonat- 
ing states at E = ±0.028(±0.025) eV for U s = 0.5(-1.0) 
eV, for instance, may not be the same, but they are of 
the same order. Therefore, we may interpret them as the 
SC pairing states influenced by the impurity. 

The response of d x 2_ y 2 wave pairing state to impu- 
rity scattering in cuprate superconductors has been in- 
vestigated intensively, and the study of the impurity ef- 
fects in superconductors has been considered as an effec- 
tive method to reveal the physics of superconductivity 2 ^. 
Fig. [5] and [5] show the LDOS in the state of d x 2_ y 2 wave 
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FIG. 4: (color online) The LDOS for anisotropic s wave pair- 
ing state for attractive scattering potential U s = —0.5 eV (a), 
U s = -1.0 eV (b), Us = -2.5 eV (c), and U s = -6.0 eV (d). 
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FIG. 6: (color online) The LDOS for d x 2_ y2 wave pairing 
state for attractive scattering potential U a = —0.5 eV (a), 
Us = -0.8 eV (b), U s = -1.0 eV (c), U a = -1.5 eV (d), 
Us = -2.0 eV (e), and U a = -6.0 eV (f). 




FIG. 5: (color online) The LDOS for d x 2_ y 2 wave pairing 
symmetry for repulsive scattering potential Us = 0.5 eV (a) 
and Us = 2.0 eV (b). 



pairing for positive and negative scattering potential, re- 
spectively. It has been found that there is no in-gap 
bound state for repulsive scattering potential cases. How- 
ever, for attractive scattering interaction, the impurity- 
induced in-gap bound states exhibit a high electron-hole 
asymmetry. The locations of the impurity states move 
towards the Fermi level from the positive energy as the 
amplitude of the scattering potential increases from 0.5 
to 1.0 (eV), and shifts to the negative energy then to 
the gap-edge as the intensity of the potential approaches 
to the unitary limit. Obviously, these are impurity- 
induced bound states stemming from the sign-change of 

the d x 2_ y 2 wave SC gap. In our calculations, no evi- analysis demonstrates that the contribution of d xy or- 
bital is larger than that of d xz and d yz orbitals which in- 
dicates the significance of the multi-orbital feature of the 
Fe-based superconductors. The distribution of the order 
parameter in real-space at unitary limit of the impurity 
scattering, as plotted in Fig. shows a color mapping 



FIG. 7: (color online) Real-space distribution of the order pa- 
rameter A(i) of anisotropic s (a) and d x 2_ y 2 (b) wave pairing 
states at unitary limit of the scattering interaction. 



dent in-gap bound states are observed when the scatter- 
ing potentials are small enough as \U S \ < 0.4 eV, which 
is physically justified that the impurity-induced bound 
states would become small and vanish as the amplitude 
of scattering interaction decreases. Decomposed LDOS 
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of the cross shaped^ suppression of the order parame- 
ters which occurs prominently along the NNN directions 
for anisotropic s wave and NN directions for d x 2 _ y i wave 
pairing state, respectively. 



Different from the SC pairing states which are influ- 
enced by the impurity in the state of anisotropic s wave 
pairing, the single impurity-induced in-gap bound state, 
for the d x 2_ y 2 wave pairing symmetry, is located near 
the Fermi energy for a moderate scattering potential 
U s = — 1 .0 eV. This is reminiscent of the cuprates where 
the emergence of the impurity-induced bound states, as 
a result of the nodal d x 2_ y 2 wave pairing state, is char- 
acterized by its low-lying feature. It has been noticed 
that the amplitude of the density of the impurity-induced 
bound states is determined by the Bogoliubov quasipar- 
ticle weights \uf a \ 2 and \v™ a | 2 , and the asymmetry of 
LDOS of the bound states in Fig. [B] demonstrates a 
strong violation of the electron-hole symmetry. There- 
fore, for FeSe-based 122-type superconductors, the node- 
less d x 2_ y 2 wave pairing state leads to a peculiar bound 
states which are characterized by a strong electron-hole 
asymmetry. Particularly, only one single robust bound 
state on impurity site has been found in the gap for 
the attractive scattering interaction. This is consistent 
with the results obtained by Zhu et al£± for repulsive 
scattering potential. In our calculation, however, it has 
been found that there is no impurity-induced bound state 
for repulsive scattering interaction in the state of d x 2_ y 2 
wave pairing, which can be interpreted as the effect of 
the strong violation of the electron-hole symmetry due 
to the peculiar Fermi surface of the FeSe-based 122-typc 
superconductors where only electron pockets exist. 



IV. SUMMARY 

In summary, using a three-orbital model, we present 
a comprehensive investigation of the effects of a single 
nonmagnetic impurity in A :r Fe2-j,Se2 (A=K, Rb, or Cs) 
superconductors by means of BdG theory. The in-gap 
bound states has been found for d x 2_ y 2 wave pairing 
state. For anisotropic s wave pairing state, the observed 
resonating states are typically SC pairing states influ- 
enced by impurity. The locations of the bound states for 
d x 2_ y 2 pairing symmetry, moving towards the Fermi en- 
ergy from the gap-edge as the amplitude of the scattering 
potential increases from 0.5 to 1.0 (cV), appear uniquely 
for attractive interaction. Due to the absence of the hole 
pocket, strong violation of the electron-hole symmetry 
occurs, which results in that no impurity-induced bound 
state emerges for repulsive scattering interaction. There- 
fore, distinguishable responses to a single nonmagnetic 
impurity, in terms of LDOS on and near the impurity 
site, would be an approach to examine the pairing sym- 
metries of FeSe-based 122-type superconductors. 



Appendix: irreducible hopping matrices of the 
three-band model 

The band structure of the three-orbital model is ob- 
tained by diagonalizing the Hamiltonian in orbital rep- 
resentation. The Z?4/j point group symmetry guarantees 
that the Bloch Hamiltonian can be written in a compact 
rotational form as the following 

H = ^2^(k)M(k)ip(k) (A.l) 

k 



Additionally, the subtle mechanism of the substitution 
of Co and Ni for Fe remains a controversial question in Fe- 
based superconductors. The effective impurity potential 
induced by various dopants which can be ascribed to the 
difference of the on-site ionic potential, together with the 
modification of the Fermi surface caused by the rigid shift 
of the band structure, has been investigated from both 
theoretical and experimental perspectiv o 37 i 43 ' 44 . The iso- 
valence of Co and Fe in Ca(Fe .944Coo.o56)2As2 com- 
pound has been pointed out and a binding energy shift 
(0.25 eV) has also been observed, which reveal the com- 
plexity of the impurity nature of Co^ 4 -. Inspired by this 
finding, one may anticipate that the doping of Ni in 
iron pnictide or chalcogenide superconductors may be 
another promising alternative since the effective impu- 
rity potential, as reported by Nakamura et al£±, is higher 
than that of Co. The impurity-induced bound states ob- 
served in our calculations in the state of d x 2_ y 2 wave 
pairing when the scattering potential is within the range 
of U s — —0.5 ~ —6.0 eV may simulate the cases of Ni 
and Co doping in FeSe-based 122-type superconductors. 



where ip(k) = [a xz (k), a^ z (k), a' xy (k)] consists of the 
orbital-component operators and 

M(k) = K 

+ K x e ik * +RK 1 R- 1 e tk « 

+ R 2 K 1 R- 2 e- ik ' +R 3 KiR- 3 e- ik * 

+ K 2 e^ k " +k ^ + RK 2 R^ 1 e l< -- k:c+k y ) 

+ R 2 K 2 R~ 2 e^- k "- k ^ + R 3 K 2 R^e l ^- k ^ 

(A.2) 

where the point group operation is carried out by rotation 

f° 1 °\ 
R = 1 (A.3) 



The irreducible hopping subsets (in unit: eV) corre- 
sponding to on-site atomic energies, hopping along x, 
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and x + y directions are For simplicity, the spin index has been dropped. 

/ 0.00 0.00 0.00 \ 
K = 0.00 0.00 0.00 
\ 0.00 0.00 0.40 J 

/ 0.05 0.00 -0.20 \ 
K x = 0.00 0.01 0.00 

\ 0.20 0.00 0.20 / 

/ 0.02 0.01 0.10 \ 
K 2 = 0.01 0.02 0.10 
\ -0.01 -0.10 0.20 / 
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